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Abstract 

We present preliminary updated measurements of the CP-violating parameters S-j^t^ and Ctttt in 
7r"'"7r~ decays. Using a sample of 227 million T(45) — > BB decays collected with the BABAR 
detector at the PEP-II asymmetric-energy e^e~ collider at SLAC, we observe 467ib33 signal decays 
and measure S^^ = -0.30 ± 0.17 (stat) ± 0.03 (syst), and C^^ = -0.09 ± 0.15 (stat) ± 0.04 (syst). 
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1 INTRODUCTION 



In the standard model, all CP-violating effects arise from a single phase in the Cabibbo-Kobayashi- 
Maskawa quark-mixing matrix [1]. In this context, neutral B decays to the CP eigenstate vr+vr" 
can exhibit mixing-induced CP violation through interference between decays with and without 
mixing, and direct CP violation through interference between the b u tree and b d 
penguin decay processes [2]. Both effects are observable in the time evolution of the asymmetry 
between B^ and B^ decays to vr+vr", where mixing-induced CP violation leads to a sine oscillation 
with amplitude S-j^t^ and direct CP violation leads to a cosine oscillation with amplitude C^^ . In the 
absence of the penguin process, C,r7r = and /Stttt = sin2Q;, with a = arg [— while 
significant tree-penguin interference leads to Ctttt 7^ and S'tttt = v^l — C^^ sin 2aefr- The difference 
between agff and a can be determined from a model-independent analysis using the isospin-related 
decays B^ 7r±7r° and B^, B^ 7r°7r° [3, 4]. 

The Belle collaboration recently reported [5] an observation of CP violation in B^ tt+tt" 
decays using a data sample of 152 million BB pairs, while our previous measurement [6] on a 
sample of 88 million BB pairs was consistent with no CP violation. In this paper we report improved 
measurements of the CP-violating parameters S-j^tt and C^jr using a data sample comprising 227 
million BB pairs collected with the BABAR detector at the PEP-II asymmetric-energy e'^e~ collider 
at SLAC. 

2 THE BABAR DETECTOR 

The BABAR detector is described in detail elsewhere [7] . The primary components used in this anal- 
ysis are a charged-particle tracking system consisting of a five-layer silicon vertex tracker (SVT) 
and a 40-layer drift chamber (DCH) surrounded by a 1.5-T solenoidal magnet, an electromag- 
netic calorimeter (EMC) comprising 6580 CsI(Tl) crystals, and a detector of internally reflected 
Cherenkov light (DIRC) providing K-tt separation over the range of laboratory momentum relevant 
for this analysis (1.5-4.5 GeV/c). 

3 ANALYSIS METHOD 

The analysis method is similar to that used in our previous measurement of these quantities [6] . We 
reconstruct a sample of neutral B mesons {Bj-gc) decaying to final states with two charged tracks, 
and examine the remaining particles in each event to determine whether the second B meson (Stag) 
decayed as a B^ or B^ (flavor tag). The CP asymmetry parameters in B^ tt+tt^ decays are 
determined with a maximum likelihood fit including information about the flavor of Stag and the 
difference At between the decay times of the Sj-cc and Stag decays. The decay rate distribution 
/+ (/_) when Srec tt+tt" and Stag = S° (S°) is given by 

e-l^*IA 

U{At) = — [1 ± 5^^sin(AmdAt) 

4r 

^ C^^cos{AmdAt)], (1) 

where r is the mean S" lifetime and Am^ is the mixing frequency due to the neutral-S-meson 
eigenstate mass difference. 

We first perform a maximum-likelihood fit that uses kinematic, event-shape, and particle- 
identification information to determine signal and background yields corresponding to the four 
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distinguishable final states (tt+tt , K~^n , K tt"*", K~^K ). The results of this fit are described in 
Ref. [8], which reports our observation of direct CP violation in S° — ^C+tt" decays. The parame- 
tjGrs S-j^'j^ cind C-j^^ are obtained from a second fit adding B-fiavor and decay-time information, where 
the yields and Kir asymmetries for signal and backgrond events are fixed to the values obtained in 
the first fit. 

3.1 Event Selection 

We reconstruct two-body neutral-B decays from pairs of oppositely-charged tracks located within 
the geometric acceptance of the DIRC and originating from a common decay point near the in- 
teraction region. We require that each track have an associated Cherenkov-angle (^c) measured 
with at least five signal photons detected in the DIRC, where the value of 9c must agree within 4(j 
with either the pion or kaon particle hypothesis. The last requirement efficiently removes events 
containing high-momentum protons. Electrons are removed based on energy-loss measurements in 
the SVT and DCH, and on a comparison of the track momentum and associated energy deposited 
in the EMC. 

Identification of pions and kaons is primarily accomplished by including 9c as a discriminating 
variable in the maximum likelihood fit. We construct probability density functions (PDFs) for 
9c from a sample of approximately 430000 D*~^ D^n^ {D^ K~7r~^) decays reconstructed in 
data, where /ir^ tracks are identified through the charge correlation with the vr^ from the D*^ 
decay. The PDFs are constructed separately for K~, tt"'", and 7r~ tracks as a function of 
momentum and polar angle using the measured and expected values of 9c, and its uncertainty. The 
average K-ir separation, defined as the difference between the expected angles for the kaon and 
pion mass hypotheses divided by the average uncertainty, varies from 12 standard deviations (cr) 
at a laboratory momentum of 1.5GeV/c, to 2a at 4.5GeV/c. 

Signal decays are identified using two kinematic variables: (1) the difference AE between the 
energy of the B candidate in the e+e" center-of-mass (CM) frame and \/i/2 and (2) the beam- 
energy substituted mass mss = \/ {s/2 + pi ■ pbY/E? — p^. Here, ^/i is the total CM energy, 
and the B momentum pB and the four- momentum of the initial state (£'j,pi) are defined in 
the laboratory frame. For signal decays, AE and tties are distributed according to Gaussian 
distributions with resolutions of 27MeV and 2.6MeV/c^, respectively. The distribution of mss 
peaks near the B mass for all four final states. To simplify the likelihood fit, we reconstruct the 
kinematics of the B candidate using the pion mass for all tracks. With this choice, B^ 7r"'"7r~ 
decays peak near AE = 0. For B decays with one or two kaons in the final state, the AE peak 
position is shifted and parameterized as a function of the kaon momentum in the laboratory frame. 
The average shifts with respect to zero are —45 MeV and —91 MeV, respectively, and this separation 
in AE provides additional discriminating power in the fit. We require 5.20 < hies < 5.29GeV/c^ 
and |A£'| < 150 MeV. The large sideband region in tties is used to determine background-shape 
parameters, while the wide range in AE allows us to separate B decays to all four final states in 
the same fit. 

We have studied potential backgrounds from higher-multiplicity B decays and find them to 
be negligible near AE = 0. The dominant source of background is the process e+e" QQ {q = 
u,d,s,c), which produces a distinctive jet-like topology. In the CM frame we define the angle 9s 
between the sphericity axis [9] of the B candidate and the sphericity axis of the remaining particles 
in the event. For background events, | cos peaks sharply near unity, while it is nearly flat for 
signal decays. We require |cos 6*51 < 0.8, which removes approximately 80% of this background. Ad- 
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Table 1: Average tagging efficiency e, average mistag fraction w, mistag fraction difference Aw = 
w{B^) — w{B^), and effective tagging efficiency Q for signal events in each tagging category. The 
quantities are measured in the .Bflav sample. 



Category e(%) w {%) Aw{%) Q {%) 

Lepton 9.6 ± 0.1 3.4 ± 0.5 0.1 ± 0.8 8.3 ± 0.3 

Kaon I 16.7 ± 0.1 8.6 ± 0.5 -1.8 ± 0.8 11.4 ± 0.3 

Kaon II 19.3 ± 0.1 20.0 ± 0.5 -3.0 ± 0.9 7.2 ± 0.2 

Inclusive 20.1 ± 0.1 30.7 ± 0.6 -4.4 ± 0.9 3.0 ± 0.2 

Untagged 34.3 ± 0.2 

Total Q 29.9 ± 0.5 



ditional background suppression is accomplished by including the Fisher discriminant J- described 
in Ref. [6] as one of the variables in the maximum likelihood fit. 

3.2 B-Flavor Identification and Decay-Time Reconstruction 

We use a multivariate technique [10] to determine the flavor of the Stag meson. Separate neural 
networks are trained to identify primary leptons, kaons, soft pions from D* decays, and high- 
momentum charged particles from B decays. Events are assigned to one of five mutually exclusive 
tagging categories based on the estimated average mistag probability and the source of the tagging 
information (Table 1). The quality of tagging is expressed in terms of the effective efficiency 
Q = Z^fcefc(l — 2wk)'^, where and are the efficiencies and mistag probabilities, respectively, 
for events tagged in category k. Table 1 summarizes the tagging performance measured in a data 
sample Sflav of fully reconstructed neutral B decays to (tt"*", p~^, a^*"). The assumption of 

equal tagging efficiencies and mistag probabilities for signal tt+tt", i^+Tr", and K'^K~ decays is 
validated in a detailed Monte Carlo simulation. Separate background efficiencies for the different 
decay modes are determined simultaneously with S'tttt and Ctttt in the fit. 

The time difference At = Az/P'yc is obtained from the known boost of the e~^e~ system 
(/37 = 0.56) and the measured distance Az along the beam {z) axis between the Bj-ec and Stag 
decay vertices. We require \At\ < 20ps and a At < 2.5ps, where a At is the error on At determined 
separately for each event. The resolution function for signal candidates is a sum of three Gaussians, 
identical to the one described in Ref. [10], with parameters determined from a fit to the -Bflav sample 
(including events in all five tagging categories). The background At distribution is modeled as the 
sum of three Gaussian functions, where the common parameters used to describe the background 
shape for all tagging categories are determined simultaneously with the CP parameters in the 
maximum likelihood fit. 

3.3 Mciximum Likelihood Fit 

We use an unbinned extended maximum likelihood fit to extract CP parameters from the -Bj-ec 
sample. The likelihood for candidate j tagged in category k is obtained by summing the product 
of event yield rii, tagging efficiency ei^k, and probability Vi^k over the eight possible signal and 
background hypotheses i (referring to tt+tt", K~^Tr~, K~n~^, and K~^K~ combinations). The 
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extended likelihood function for category k is 



exp 



E 



n 



(2) 



The yields for the Kir final state are parameterized as npc±^T = fiK-w (1 T -A-k-k) /2, where Ak-k 
is the direct-CP-violating asymmetry [8]. The probabilities Vi,k arc evaluated as the product of 
PDFs for each of the independent variables Xj = {m-Es, ^E, J^,9^ ,9^ , At}, where 9^ and 9~ are 
the Cherenkov angles for the positively- and negatively-charged tracks. The At PDF for signal 
tt'^tt" decays is given by Eq. 1, modified to include the Wk and Awk for each tag category, and 
convolved with the signal resolution function. The At PDF for signal Ktt decays takes into account 
mixing and the correlation between the charge of the kaon and the flavor of -Stag- 
There are 46 free parameters in the fit: 

• 12 parameters describing the background PDFs for niES, ^E, and J^; 

• 8 parameters describing the background At PDF; 

• 12 background flavor-tagging efficiencies; 

• 12 background flavor-tagging efficiency asymmetries; and 

• '^TTTT cind Ctttt . 

The signal and background yields and Kir asymmetries were determined in a separate fit that 
does not use flavor-tagging or At information [8]. Out of a fitted sample of 68030 events, we find 
iT'TTTv = 467 ± 33, uktt = 1606 ± 51, and ukk = 3 ± 12 decays, and measure Aktt = —0.133 ± 0.030, 
where all errors are statistical only. We fix r and Am^ to their world-average values [11]. The 
total likelihood jC is the product of likelihoods for each tagging category, and the free parameters 
are determined by maximizing the quantity ln>C. 



4 PHYSICS RESULTS 

The fit to the Bj-ec sample yields 

S^^ = -0.30 ± 0.17 (stat) ±0.03 (syst), 
C^^ = -0.09 ±0.15 (stat) ±0.04 (syst). 

The correlation between St^t^ and C\.„ is —1.6%, and the correlations with all other free param- 
eters are less than 1%. These preliminary results are consistent with our previously published 
measurements [6], and are combined with our measurements of the branching fractions for the 
isospin-related decay modes tt^tt^ and tt^tt^ to determine model-independent bounds on a (see 
Ref. [4]). 

Figure 1 shows distributions of ttt-es and AE for events enhanced in signal tt+tt" decays. We 
apply additional requirements on probability ratios based on all PDFs except the variable being 
plotted. The solid curves are projections of the maximum likelihood fit for the sum of signal and 
all background components, while the dashed curves indicate the sum of qq and B^ K~^tt~ 
cross-feed background. Figure 2 shows distributions of At for events with i?tag tagged as B^ or 
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Figure 1: Distributions of (left) mss and (right) AE for events (points with error bars) enhanced 
in signal tt^tt" decays using additional requirements on probability ratios. Solid curves represent 
projections of the maximum likelihood fit, dashed curves represent qq and Ktt — > vrvr cross- feed 
background. 




-6 -4 -2 2 4 



At (ps) 



Figure 2: Distribution of the decay-time difference At for events enhanced in signal B ■k'^tt~ de- 
cays using additional requirements on probability ratios. The top two plots show events where -Btag 
is identified as (a) (ra^o) or (b) B^ {n^o)- Solid curves indicate the projection of the maximum 
likelihood fit including signal and background, while the dashed curves show the contribution from 
background events, (c) The asymmetry (points with errors), defined as (n^o — n^o) / (n^o + n^o), 
for different bins in At, and the projection of the full fit (solid curve). 
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B^, and the asymmetry as a function of At for signal tttt decays selected with requirements on 
probability ratios including all PDFs except At. 

As a consistency check on the At resolution function, we perform a B^-B^ mixing study using 
the large number of Ktt signal decays in the Srcc sample. Floating r and Arrid along with Stttt, 
Ctttt) and Ak-k-, wc find values consistent with the world averages (r = 1.60 =b 0.04 ps and Arud = 
0.523 ± 0.028 ps^^), and the CP parameters are consistent with the nominal fit results. 

5 SYSTEMATIC STUDIES 

Table 2 summarizes the contributions to the total systematic uncertainty from the dominant 
sources. These include imperfect knowledge of the PDF shape parameters; the B-fiavor-tagging pa- 
rameters (Table 1); the alignment of the SVT; the event-by-event beam-spot position; the potential 
effect of doubly Cabibbo-suppressed decays of the i?tag meson [12], and the B lifetime and mixing 
frequency. In addition, to confirm that we are sensitive to non-zero values of S'tttt and Ctttt, we 
fit a large sample of Monte-Carlo simulated signal decays with large values of the CP parameters. 
The fit results are consistent with the generated values, and we assign the sum in quadrature of 
the statistical uncertainty and the difference between the fitted and generated values as a conser- 
vative systematic error accounting for potential bias. The effect of uncertainty on the signal and 
background yields and Kt: asymmetries is negligible for both 5^,^ and C^tt- The total systematic 
uncertainty is calculated by summing in quadrature the individual contributions. 

Table 2: Summary of systematic uncertainties on 6*7^77 and Ctttt (sec text for details). The total 
uncertainty is calculated as the sum in quadrature of the individual contributions. 



Source 






PDF parameters 


0.017 


0.018 


i?-flavor identification 


0.005 


0.015 


SVT alignment 


0.010 


0.002 


Beam spot 


0.010 


0.010 


Tag-side interference 


0.008 


0.023 


T50 and Arrid 


0.001 


0.004 


Potential bias 


0.013 


0.007 


Total 


0.027 


0.035 



6 SUMMARY 

In summary, we present preliminary updated measurements of the CP-violating asymmetries St^t^ 
occuring in the time distributions of B^ tt+tt decays. We find S'tttt = — 0.30ib0.17ib0.03 
and Ctttt = —0.09 it 0.15 ± 0.04, which are consistent with our previous measurements and with 
the hypothesis S^^t^ = Ctttt = 0. These results do not confirm the observation of large CP violation 
reported in Ref. [5]. 
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